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Abstract. We study ttie linear dynamic mechanical response of monodomain nematic liquid crystalline 
elastomers under shear in the geometry that allows the director rotation. The aspects of time-temperature 
superposition are discussed at some length and Master Curves are obtained between the glassy state and the 
nematic transition temperature T„i. However, the time-temperature superposition did not work through 
the clearing point Tni, due to change from the "soft-elasticity" nematic regime to the ordinary isotropic 
rubber response. We focus on the low-frequency region of the Master Curves and establish the power-law 
dependence of the modulus G' oc uj°' . This law agrees very well with the results of static stress relaxation, 
where each relaxation curve obeys the analogous power law G' oc t~°' in the corresponding region of long 
times and temperatures. 

PACS. 83.80.Dr Elastomeric polymers - 61.30.-v Liquid crystals. - 83.10.Nn Polymer dynamics 



1 introduction 

For the last two decades, thermotropic liquid crystalline 
polymers (LCP) have been a subject of increasing activity 
due to their fascinating potential in technical applications, 
as well as some fundamentally different properties in com- 
parison to ordinary liquid crystals. Typical LCP consist 
of a polymer backbone and mesogenic groups. The back- 
bone is usually made of siloxane or acrylic chains and the 
mesogenic groups are either grafted (end-on or side-on) 
using a flexible spacer, thus making a side-chain LCP, or 
incorporated directly into the polymer backbone (main- 
chain LCP). When the polymer is chemically crosslinked, 
it forms a three-dimensional permanent network connected 
with mesogenic moieties in what is called Liquid Crys- 
talline Elastomers (LCE) [Q, or gels if the network is 
swollen by a solvent. Many new physical results have been 
reported for LCE, but perhaps the most important nov- 
elty, unique for nematic LCE materials, is the effect of 
soft elasticity, where imposed strains could be completely 
accommodated without any (or with very low) elastic re- 
sponse due to the relaxation of internal microstructure. 
These unusual mechanical properties of LCE can be di- 
rectly applied in technology, using LCE as anomalous damp- 
ing materials, as well as artificial muscles or mechanical 
transducers. 

The early research of LCP rheology, focused mostly on 
the viscosity of polymers, was summarized in Many 
interesting results have come to light once the mesogenic 
groups were incorporated into linear or branched poly- 
mer systems. The classical Leslie-Ericksen theory of low- 
molecular-weight liquid crystals is, in most of the cases, 
not applicable to LCP systems because of the underly- 
ing polymer dynamics making the response distinctly non- 



Newtonian. Hanus et al. [g| measured dynamic-mechanical 
properties of a weakly crosslinked main-chain LCP hav- 
ing crystalline, smectic B, smectic A and isotropic phases. 
They only observed a clear step in the linear shear mod- 
ulus G" at the transition temperature from smectic A to 
the isotropic state. In the study of LCP having smectic 
C, smectic A, nematic and isotropic phases Q, Pakula 
reported that decreases in the linear modulus were de- 
tected at the smectic/nematic transition as well as the 
glass/smectic transition. Also interestingly, although it 
was not discussed theoretically at that time, the increase 
of the modulus could be seen during the transition from 
the nematic to the isotropic phase. The issues of time- 
temperature superposition in LCP rheology have been ac- 
tively discussed by several groups such as Colby et al. 
1), Gallani et al. @,|, Four maux-Demange et al. ^ and 
Weilepp et al. All of these works also reported the 
characteristic power-law frequency dependence of the stor- 
age modulus G'{uj). Colby et al. used a sample possess- 
ing three phases (isotropic, nematic and smectic A) and 
concluded that the viscoelastic response is very sensitive 
to smectic-nematic and smectic-isotropic transitions, but 
insensitive to the nematic-isotropic transition. The time- 
temperature superposition worked only at the nematic- 
isotropic transition. Approaching the glass transition, G'{u>) 
obeys a power law of G' cx w^-''^. The authors suggested 
that this power law might be connected to the large poly- 
dispersity of their sample H, which was also argued by 
Fourmaux-Demange et al. |^ for their different LCP sys- 
tem. Gallani et al. measured dynamic-mechanical response 
of several LCP materials and also found a power law be- 
haviour of the storage modulus G", with the exponent of 
0.44, 0.6, 0.67 or 0.75 depending on the sample. Here 
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again, the time-temperature superposition worked only 
during the nematic-isotropic transitions. This was attributed 
|^,|| to the transient elastic clusters of macroscopic size 
(a few tens of micrometers), floating in a viscous medium. 
One may recall that a classical Rouse dynamics of polymer 
chains would lead to a square-root dependence G' oc w"-^ 

In this paper we shall find that the power law depen- 
dence is also effective in crosslinked LCE materials. In 
their parallel study of LCE, Gallani et al. and Weilepp et 
al. H,]^ also clearly stated that the rheological response 
did not change significantly when crossing an isotropic- 
nematic phase boundary, whereas there has been a crucial 
change at the nematic-smectic A transition. The power- 
law frequency dependence of the elastic moduli has again 
been reported, with G' oc 0;°'^ in the smectic A phase. In 
the isotropic phase, the exponent was found to be 0.5, sim- 
ilar to the Rouse model prediction. A more detailed recent 
study of nematic LCE rheology by Stein et al. also re- 
ported G' oc Ld^-^ in a polydomain nematic elastomers, but 
could not obtain Master Curves for monodomain, aligned 
nematic materials. Working with such monodomain ne- 
matic rubbers, Clarke et al. reported a significant 
difference in response between polydomain samples (where 
the nematic director is disordered on the scale above few 
microns) and monodomain nematic rubbers (with the di- 
rector is macroscopically aligned along the shear direction, 
or perpendicular to it). A substantial decrease in the low- 
frequency storage modulus G' was registered as the tem- 
perature was decreasing through nematic-isotropic transi- 
tion. This unusual effect has been theoretically described 
by a continuum theory of linear viscoelastic response in 
oriented monodomain nematic elastomers jl^ . 

Continuous stress relaxation is a related subject, which 
offers a different angle of approach to the complex prob- 
lem of polymer dynamics exploring the time, rather than 
the frequency domain. Chasset and Thirion, in their clas- 
sical paper [nSl, described a power law decay in ordinary 
crosslinked polymer networks such as the isoprene rubber 
(IR) and the butadiene-styrene copolymer rubber (SBR). 
They found that, after application of a constant strain, 
the stress response decreases with time, reaching at long 
times a characteristic power-law regime oc 1/^0.1-0.15 
very long times. In the linear regime (at small enough 
deformations), in an incompressible material (when the 
extension modulus is simply three times the shear modu- 
lus G), this translates directly into the relaxation of the 
effective modulus G{t). Many subsequent experiments on 
stress relaxation, performed by different groups, have sup- 
ported this long-time power law decay with a small expo- 
nent. Theoretical analysis of applying the reptation 
ideas for a network with long free dangling chain ends, has 
shown that the retraction of these chains could account for 
the observed stress relaxation law a oc with a small 
exponent a ^ 0.1 — 0.15. 

In liquid crystalline networks, the early results shown 
a different response. The problem in such experiments is 
always the need to access very long times, when the main 
transient processes in the polymer network have relaxed 



and only the slowest mode remains. A peculiar dual regime 
of long-time relaxation in a polydomain siloxane nematic 
LCE has been reported by Clarke et al. . A character- 
istic time i*^3000s separated a regime of fast power-law 
relaxation, G oc t~^'^ or higher, and a very slow relaxation 
at the later stage, fitted in |]l^ with an inverse-logarithmic 
law (T~[l -I- ln(i/i*)]~^. An even faster power-law decay 
was reported by Hotta et al. |l^ in their study of acrylic 
nematic LCE, G oc t^"-^^ in the early stage of stress re- 
laxation. At very long times, they find the different slow 
power law decay with the exponent ^0.15, similar to the 
one observed in natural crosslinked rubber by Chasset and 
Thirion. 

From this brief summary, it is clear that the power- 
law dependence of linear response function on time or fre- 
quency is a recurrent feature of LCE rheology. On the 
other hand, it is equally clear that the response is non- 
universal and the particular features of the response de- 
pend on the material composition and also on the ne- 
matic director texture and orientation with respect to 
the applied strain. In this paper we describe a simulta- 
neous study of dynamic-mechanical and stress-relaxation 
responses for the two well-characterised nematic LCE ma- 
terials. One sample is a "classical" siloxane side-chain ne- 
matic rubber first prepared and described in great detail 
by Finkelmann et al. and others ||l|,|o[[l^,|l|l. The other 
material had a high content of main-chain nematic poly- 
mer strands, which makes the effective chain anisotropy 
very high but also leads to a very slow dynamics. We 
shall examine in some detail the issues of time-temperature 
superposition, which is necessary to extrapolate the re- 
sponse function into the regions of time or frequency not 
practically accessible on experiment. We find that the re- 
sults of two such different experimental techniques are, in 
fact, in a good agreement: the observed frequency power 
laws G' oc uj"" are matched by the corresponding relax- 
ation laws G oc t"", if the corresponding ranges of time 
and frequency are examined. 



2 Material Characterisation 

Two nematic elastomers were prepared in our group fol- 
lowing the synthetic procedure pioneered by Finkelmann 
et al. The polymer backbone was a standard poly- 
methylhydrosiloxane with approximately 60 Si-H units 
per chain, obtained from ACROS Chemicals. The pendant 
mesogenic groups were 4'-methoxy phenyl-4-(l-buteneoxy) 
benzoate (MBB), as illustrated in Fig. |l|, attached to the 
backbone end-on via the hydrosilation reaction in the pres- 
ence of a commercial platinum catalyst COD, obtained 
from Wacker Chemie. The first monodomain nematic elas- 
tomer was prepared following the two-step crosslinking re- 
action |T^, with 1,4 di(ll-undecene) benzene (HUB), a 
small fiexible difunctional molecule deemed to have rela- 
tively minor effect on the overall mesogenic properties of 
the liquid crystalline polymer (synthesized in the house). 

The second network was prepared following the similar 
procedure, but using di-vinyl terminated polymer chains 
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Mesogenic monomer 



(MBB) 



o ^ 10% J ^[ crosslmkers: (i) or (ii) 



(a) 



(i):11UB 



(b) 



(ii) : MC 



Fig. 1. Schematic illustration of the materials used in this 
work, (a) Siloxane backbone chain with Si-H groups reacting 
with 90 mol% mesogenic phenyl-benzoate side groups, MBB 
and 10 mol% of di- vinyl crosslinking groups: (b) flexible small- 
molecule 1,4 alkeneoxybenzene, HUB (resulting in the SiH ma- 
terial), and the main-chain nematic polymer of l-biphenyl-2- 
phenyl butane, MC (giving the SiMC material). 



of a-{4-[l-(4'-{ll-undecenyloxy} biphenyl)-2-phenyl] butyl) 
a;-(ll- undecenyloxy) poly-[l-(4- oxydecamethyleneoxy)- 
biphenyl -2-phenyl] butyl (MC) that themselves make a 
mam- chain LCP [E2 . In our case, the crosslinking chains 
had ~75 rod-hke monomers between the terminal vinyl 
groups (determined by GPC, polydispersity '^2). As a re- 
sult, the properties of a nematic LCE are dominated by 
these long main-chain LCP strands connecting the silox- 
ane side-chain molecules (as the Table 1 indicates, the 
overall MC content is ~58%). In both cases, the crosslink- 
ing density was kept at 10 mol% of the reacting bonds 
in the siloxane backbone, so that on average, each silox- 
ane chain has nine pendant mesogenic groups between 
crosslinking sites. 

Both elastomers were prepared from the same batch 
of backbone and side-group mesogenic materials, with the 
same relative concentration of crosslinking groups. 

Table 1 summarizes chemical and thermal character- 
istics of both materials. Equilibrium transition tempera- 
tures were determined on a Perkin Elmer Pyris 1 differ- 
ential scanning calorimeter (DSC), extrapolating to low 
cooling rates, and the nematic phase identified by polariz- 
ing optical microscopy and X-ray scattering. The average 
degree of chain anisotropy was identified by observing the 
uniaxial thermal expansion of LCE samples as a function 
of temperature. Thermal expansion measurements were 
made by suspending the samples, without load, in a glass- 
fronted oven and measuring the variation in natural length 
of the samples with temperature, L(T); see 21 for further 
details. 



Sample 


HUB 


MC 


SC wt 


r 




J- m 


SiH 


10% 


0% 


87% 


1.5 


3°C 


86° C 


SiMC 


9% 


1% 


42% 


2.8 


2°C 


107° C 



Table 1. Proportions (in mol%) of crosslinkers HUB and MC 
in the overall crosslinking composition (of the fixed total of 
10%), the corresponding volume fraction of side-chain meso- 
genic groups (in wt%), the average chain anisotropy r at room 
temperature, temperatures of glass and nematic-isotropic tran- 
sitions. The glass transition temperatures are approximate, 
with an error of at least ±5°C. 





Sandwich 


Tape 




Thickness 


Area 


Thickness 


Area 


SiH 
SiMC 


0.25 mm 
0.15 mm 


79 mm^ 
(tt X 5^) 
36 mm^ 
(tt X 3.4^) 


2 mm 
2 mm 


1.9 mm^ 
(7.4 X 0.25) 

1 mm^ 
(6.7 X 0.15) 



Table 2. Samples and their dimensions for the simple shear 
experiment in two mounting geometries, see text. 



3 Dynamic Mechanical Analysis 



In this paper, we used the dynamic mechanical appara- 
tus (DMA) Viscoanaliser VA4000, from Metravib RDS. In 
the previous work |13,^ we have already examined the 
anisotropy of the shear response depending on the orien- 
tation of the nematic director, see the sketch in Fig. ||(a). 
Only the rotation-inducing geometry reflects the effects of 
soft elasticity and here we concentrate on this particular 
case. The way in which the samples were mounted in the 
DMA is shown in Fig. |(b). It was important to choose 
the correct mounting for a given temperature range be- 
cause of the great variation in the magnitude of response 
between the low-temperature (or high-frequency) glassy 
regime and the high-temperature rubber plateau. We used 
two different ways of sample mounting: The more stiff 
"sandwich" geometry, with each sample having the large 
area of contact with the clamps and the small thickness, 
was used to collect weaker signals from softened samples 
at higher temperature (in our DMA setup the thin-film 
sample had a round shape). The "tape" mounting geom- 
etry, in which the rectangularly cut thin-film sample is 
clamped sideways with only the small cross-section area, 
results in smaller forces and allows to measure the mod- 
ulus of very rigid samples up to the glass plateau. The 
Table 2 gives the dimensions of typical sample mounting. 

A geometry of uniaxial extension, more commonly used 
in studies of equilibrium stress-strain relaxation in rub- 
bers, is less appropriate for an oscillating regime because 
of slow relaxation and incomplete sample recovery on each 
cycle. In the chosen simple shear setup, the middle part of 
the sample holder is forced to move in an oscillating fash- 
ion, maintaining the small shear strain - between 0.003 
and 0.004 (that is, 0.3-0.4%). This value was chosen from 
a separate strain-ramping experiment, finding the range 
of strains that still return a linear response in the whole 
of the chosen frequency range. 
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load cell 




(a) (b) 

Fig. 2. (a) The geometry of simple shear with two principal 
orientations of the initial director n, labelled D (for n along 
the shear displacement) and V (vorticity). No director rota- 
tion, hence no soft elasticity, occurs in the V-geometry jl^. 
(b) The symmetric sample mounting in the DMA device, in 
the Z)-geometry. See text for the sample dimensions used in 
different setups. 

3.1 Temperature Scans 

The materials were studied over a wide temperature range 
encompassing the glass and the nematic/isotropic transi- 
tions. Typically the temperature scans ranged from — 50°C 
to 130°C, at a number of fixed frequency values. Fig- 
ure ^ shows the results for the storage modulus G' at 
just two frequencies, 0.01 and 100 Hz, for the SiH sam- 
ple. The storage modulus G' , which is the real part of 
the linear complex modulus G* {ui) and the loss factor 
tan S = G" /G' (with G" the imaginary part of G*{uj) - 
the loss modulus) were plotted against the temperature. 
As we mentioned above, it is technically difficult to obtain 
a clean mechanical signal from the sample that undergoes 
a change in stiffness of up to 10^ — > lO^Pa. Therefore, we 
used the two ways of mounting the samples, the "sand- 
wich" at high and the "tape" at low temperatures. One 
can see in Fig. ^, as well as in Fig. ^ for the SiMC material, 
that the two portions of the full curve, for the tape and 
sandwich modes, are connecting quite successfully with- 
out any additional adjustment (the data were discarded if 
the vertical error was reaching 15%, on the assumption of 
a mounting fault). This indicates good reliability of DMA 
measurements for both the storage modulus G' and loss 
factor tan 5. 

The glass plateau modulus for both SiH and SiMC 
materials reaches up to nearly 10^ Pa, which is a reason- 
able value in polymeric materials. Between 0°C and 10°C, 
both SiH and SiMC undergo the glass/nematic transition, 
showing a dramatic decrease in G" on cooling, and the as- 
sociated steep rise in tan 5. The storage modulus decreases 
down to ~ 10^ Pa which is an expected order of magnitude 
in rubbery phases. However, both materials show a sub- 
stantial drop in G' in the nematic phase, below Tni- This 
effect, most pronounced at low frequencies, is the dynamic 
soft elasticity: the reduction in the modulus and the rise in 
internal mechanical dissipation caused by the underlying 
director rotations induced by shear The relative 

depth of the drop in G' is higher in SiMC, the material 
with a higher chain anisotropy due to its main-chain LCP 
component. Above T„i, the storage modulus returns to its 
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Temperature ( C) Temperature ( C) 

(a) (b) 

Fig. 3. Storage modulus G' (circles) and loss factor tan 5 
(diamonds) as functions of temperature at frequencies of 0.01 
(a) and 100 Hz (b). The graphs combine the data obtained 
from the tape (open symbols) and the sandwich (filled symbols) 
sample mounting geometry. Results for SiH. 




-50 50 100 -50 50 100 



Temperature ( C) Temperature ( C) 

(a) ft)) 

Fig. 4. Storage modulus G' (circles) and loss factor tan 5 
(diamonds) as functions of temperature at frequencies of 0.01 
(a) and 100 Hz (b). The graphs combine the data obtained 
from the tape (open symbols) and the sandwich (filled symbols) 
sample mounting geometry. Results for SiMC. 

(higher) value characteristic of ordinary isotropic rubbers. 
At very low frequencies G' is slightly rising with temper- 
ature, revealing that the sample is in a crosslinked rub- 
bery state with its equilibrium shear modulus G ~ n^kBT 
(where Ux is the crosslink density and fcs is Boltzmann 
factor). 

3.2 Frequency Scans 

Figures || and || show the same data for the storage mod- 
ulus G", but this time plotted against frequency at a set 
of constant values of temperature. In our experiments, the 
frequency was changing between 0.01 to 900 Hz, measured 
45 points in total. However, we were obliged to discard 
higher frequencies: the DMA data sets must pass a series 
of tests on internal consistency and failing these tests in- 
dicates one of several possible problems with acquisition 
(most common is the internal resonance of the mechanical 
frame) . Only the reliable range of results is presented in 
Figs. ^ and 

The range of storage modulus variation with frequency 
matches well with the temperature scans, from 10^ Pa 
down to 10"* — 10^ Pa. As in the temperature scans, the 
average value of G' is rising on heating, when each sam- 
ple approaches its isotropic rubbery phase and the meso- 
genic effects weaken. In each of the figures, the plot (b) 
shows the expanded version of the high-temperature re- 
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Frequency (Hz) 



Fig. 5. Storage modulus G' as a function of frequency at tem- 
peratures ranging from — 10°C to 130°C, labelled on plots; 
results for SiH. Plot (a) shows the set of lower temperatures, 
including the glass and the nematic rubber plateau regimes and 
clearly suggesting the time-temperature superposition. Plot (b) 
shows the higher temperatures, below and above Tni, indicat- 
ing the rise in the rubber plateau value in the isotropic phase: 
the dynamic signature of soft elasticity. 




frequency (Hz) 



FreqLiencv (Hz) 



(a) 



(h) 



Fig. 6. Storage modulus G' as a function of frequency at tem- 
peratures ranging from — 5°C to 130°C, labelled on plots; re- 
sults for SiMC. Plot (a) shows the set of lower temperatures, 
including the glass and the nematic rubber plateau regimes and 
again suggesting the time-temperature superposition. Plot (b) 
shows the higher temperatures, below and above Tni, same as 
in Fig. 1(b). 



suits, clearly indicating the effect of dynamic soft elastic- 
ity. In addition, one can see for the SiMC material that the 
low-frequency rubber plateau is not achieved even at high- 
est temperatures: the modulus continues to decrease with 
decreasing frequency indicating the on-going mechanical 
relaxation processes in the sheared rubbery network. It 
is in contrast with the SiH results, where one can clearly 
identify the frequency-independent rubber plateau in both 
the nematic and the isotropic phases. 

Below we shall discuss this extremely slow relaxation 
of the MC-containing elastomers and match these results 
with the findings of static stress relaxation. We attribute 
this slowing down to the hairpin folds of the backbone of 
rods with flexible spacers, which has been demonstrated 
p5| to freeze the dynamics even of non-crosslinked poly- 
mer melts. 



3.3 Time- Temperature Superposition 

The frequency scans presented in Figures 5|(a) and ^(a) in- 
vite the time-temperature superposition 1 24 . At the same 
time, one immediately discovers from the corresponding 
(b)-plots that such superposition will not be possible as 
the monodomain nematic rubber sheared in the "soft" 
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Fig. 7. Time-temperature superposition of frequency scans for 
SiH, at different temperatures, (a) The Master Curve built at 
Trof = 22° superposing the data at temperatures labelled 
on plot. Annotations also show the expected Rouse time and 
the power-law exponent of the approach to the glass transition, 
(b) Logarithm of shift factors plotted against temperature. The 
curve was fitted with the WLF eq. at lower temperatures. 
Square symbols refer to the discussion of Section 4 below. 




1|TCC)| 



'0.07 
,Q-ie ,(j-io .,0-4 
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I i, i 

10= 10' 



280 300 320 340 360 

Abs Temperature (K) 



Fig. 8. Time-temperature superposition of frequency scans for 
SiMC, at different temperatures, (a) The Master Curve built 
at Trof = 20°C, superposing the data at temperatures labelled 
on plot. Annotations also show the expected Rouse time and 
the power-law exponents. The plateau level is not reached even 
at such low extrapolated frequencies, (b) Logarithm of shift 
factors plotted against temperature. The curve was fitted with 
the WLF equation at lower temperatures. Square symbols refer 
to the discussion of Section 4 below. 



£>-geometry undergoes the nematic transition. We assert 
that this conclusion can only be made for materials which 
show a clean and pronounced soft decrease in G": in a 
frustrated polydomain elastomer, or in an overcrosslinked 
network, the required nematic director rotation may be 
impeded by internal constraints and, therefore, the char- 
acteristic nematic soft elasticity masked. In that case, one 
might expect the ordinary time-temperature superposi- 
tion to be valid, same as in ordinary isotropic polymer 
systems. 

Time-temperature (or, in this case, frequency-temperature) 
superposition is a very interesting, practically useful, but 
also theoretically ill-defined approach to data analysis in 
dynamic mechanical measurements. The empirical algo- 
rithm proposed by Williams, Landel and Ferry (WLF) 
[ p4[ allows extrapolating the results into the time or fre- 
quency ranges far beyond experimentally accessible. As 
traditional in this field, we attempted such a superposi- 
tion of our data for nematic elastomers. The storage mod- 
uli G taken in the series of frequency scans at different 
temperatures are shifted along the frequency axis until 
each data sets superimposes with the previous (with no 
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additional vertical shifting). The multiplicative shift fac- 
tor, Lo' = axi-o, is strongly temperature dependent. An- 
other important factor is the choice of reference temper- 
ature Tref at which the real frequency values are taken. 
WLF suggested using the glass transition temperature as 
Trof, however, we decided against this for the following 
reasons. First of all, the presented results show that the 
dynamic glass transition is, in fact, difficult to pinpoint as 
the transformation occurs over a range of temperatures. 
Secondly, no additional information is extracted from the 
choice of T^d anyway; all it affects is the overall shift of 
the Master Curve along the frequency axis. Our choice of 
Trof as the ambient room temperature was dictated by the 
possibility of determining the key parameter of polymer 
chains, the Rouse frequency = ^ at the relevant Trcf , 
from the position of the dynamic glass transition on the 
frequency axis. In both materials, we find the Rouse time 
tr ~ 10~^ — 10~^s (which of course must be regarded as 
a very crude estimate). 

Figures |7|(a) and ^(a) show the results of superposi- 
tion for the SiH and SiMC samples, respectively, in the 
form of the Master Curves. The procedure is only applica- 
ble at the temperatures below the nematic transition T„i. 
Each graph also shows one data set G'{uj), appropriately 
shifted, at a high temperature, in the isotropic phase. The 
values are much higher than on the low-frequency nematic 
plateau and, of course, no shifting along the w-axis would 
superimpose them. 

The very existence of Master Curves indicates that 
there is a unique physical process of mechanical relaxation, 
which spans the range of times, or frequencies that is far 
beyond experimentally accessible. Our measurements, at 
each temperature, are only able to detect a small portion 
of the whole process, but the time-temperature superpo- 
sition suggests a way of building the whole picture from 
small pieces. An important part of this analysis is the tem- 
perature dependence of shift factors ut , which is plotted 
in Figs. 0(b) and ||(b). The empirical WLF relationship is 
based on the ideas of glass transition and has been suc- 
cessful in describing the properties of simple homopolymer 
melts: 

Ci(T-rrof) , , 

^ ' c. + (T-r..eO ' 

where Ci and C2 are the constants, called the WLF co- 
efficients. (In more complex polymers, and in crosslinked 
networks, deviations from this equation have been sub- 
stantial). In our case. Figures |^(b) and ||(b), the shift fac- 
tors follow the WLF relation very well at lower tempera- 
tures (below Trof), with the WLF coefficients Ci w 17.4 
and C2 ~ 102 for both materials,]^ very close to the val- 
ues found in ordinary isotropic polymer systems. This is 
perhaps not surprising when one recalls that the empirical 
WLF relationship is "designed" to describe the dynamic 
glass transition region. 

However, the shift factors deviate from the classical 
WLF behavior very strongly as soon as the higher temper- 
atures (and, correspondingly, long relaxation times) are 

^ Recall that both SiH and SiMC have been analysed with 
nearly the same reference temperature, around 20°C. 



considered. Here, the two materials behave quite differ- 
ently: SiH only allows the superposition of frequency scans 
with much smaller shifts ot, while the Master Curve for 
SiMC demands much greater shifts. One should note that, 
in practice, the procedure of shifting of plots along the fre- 
quency axis to achieve superposition is very sensitive to 
the value of each shift factor ar and, accordingly, the er- 
rors in Figs. |7|(b) and ||(b) are very small. It is more likely 
that the noise in the ut data is due to uncertainties in 
DMA temperature or modulus measurements. 

Finally, we note that the region of the dynamic glass 
transition, in both Master Curves for SiH and SiMC, fol- 
lows a clear power law G" cx uj'^. The plots indicate the 
non-universal values of the exponent, similar to many pre- 
vious observations. The very low frequency (long-time) 
tail of continuing relaxation in SiMC may appear to fol- 
low another power law with a much smaller exponent of 
0.07 also labelled in Fig. ||(a). However, we are inclined to 
believe that this is an ambiguity of the complex transient 
regime: in the previous studies of static stress relaxation 
[0|l^ we were distinctly unable to distinguish between a 
very low power-law exponent and an even slower logarith- 
mic relaxation regime. Since the material never reached 
the equilibrium, whether in a real experiment or in extrap- 
olation, the identification of the longest relaxation regime 
remained impossible. 



4 Static Stress Relaxation 

Stress relaxation curves were measured at several differ- 
ent temperatures at a constant extensional strain of 0.09 
for both SiH and SiMC, following the simple technique 
described in |l^,|l^. This strain is higher than the shear 
strain used in DMA experiments (for technical reasons we 
simply could not impose a smaller extension) and the ques- 
tion of whether the response is still linear may be raised. 
We have to stress two points. First of all, the strain of 
^9% is, in fact, on the border of the initial linear regime 
of stress-strain curves measured for both SiH and SiMC, 
cf. II^. This, however, has to be regarded with caution be- 
cause it is impossible to obtain a truly equilibrium stress- 
strain relation (a certain rate of strain is always involved) 
and any extrapolation to zero strain rate would involve a 
theoretical model, which may or may not be ultimately 
successful. However, as a second point, we have observed 
that the key conclusion of this section, about the power 
laws of stress relaxation, remains valid even for higher 
strains - on and above the soft plateau when we are cer- 
tain that the stress-strain is non-linear. Therefore, in prac- 
tice, the choice of the strain step for the static relaxation 
experiments was not very important. 

For the constant-strain extension setup, each sample 
was cut with a razor blade into a rectangular shape of 
approximately 6 x 2.5 x 0.2 mm, measured by a microm- 
eter (narrow elongated strips minimize the edge effects 
during the extensional deformation). The initial sample 
cross-section area of 0.5 mm^ was used to calculate the 
nominal stress a, in units of (mN/nim^ = kPa), from 
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Fig. 9. Typical stress relaxation curves for the SiH (a) and 
SiMC (b) samples, after a step-strain of 9% is applied perpen- 
dicular to the nematic director (log-time axis). In each plot, 
the data sets for different temperatures are labelled. 



the force measured by the temperature-compensated dy- 
namometer (UF1/AD20 from Pioden Controls Ltd). The 
values of force, obtained in arbitrary units, were calibrated 
with weights at different experimental temperatures of 
30, 50, 80, 105°C. The samples were then extended with 
the engineering strain e = (L — Lq)/ Lq, increasing from 
to 0.09 within a second. Both monodomain samples, 
SiH and SiMC, were stretched perpendicular to their uni- 
form director orientation so that a soft-deformation geom- 
etry, comparable to the U-shear geometry of DMA exper- 
iments, is achieved. In each case the stress relaxation was 
measured for two days (~ lO^s). Figure ^ shows typical 
stress relaxation curves taken in the series of these mea- 
surements. The data were then fitted with several possible 
stress relaxation equations, such as the power-law, the in- 
verse logarithmic law and the exponentials. 

In Fig. ^(a), for the "faster" SiH material, one can 
note that the stress relaxation was substantial at a lower 
temperature of 31°C and is almost unnoticeable at higher 
temperatures (in fact, the data set for 82°C is not shown 
because it is hard to distinguish from the 52°C case). The 
relaxation in SiMC continues even at high temperatures, 
but its extent is clearly smaller. 

Instead of presenting and discussing the individual re- 
laxation fits and their relations, we shall follow a different 
route of analysis. If we assume that the Master Curves, 
shown in Figs, ^(a) and ^(a), are genuine and represent 
the complete relaxation dynamics in our complex poly- 
meric systems, then we should expect that the same laws 
must govern the static stress relaxation. Let us plot such a 
Master Curve against the inverse frequency, which would 
be a measure of real time. Such plots. Figs. [lO|(a) and pl](a) 
look reasonable, as one would expect the relaxing modu- 
lus to behave in real time, in double-logarithmic scale. Of 
course, one cannot assign anything more than a qualitative 
meaning to this representation because, let us not forget, 
the G*{uj) and G{t — t') are related by the Fourier trans- 
formation and merely mapping G(uj) onto the inverse fre- 
quency is not a rigorous operation. It does, however work 
for the power-law dependence: from the basic definitions 

a(t) = £ G{t^t')-^e{t')dt' ■ (2) 
cr(ti>) iLo G* {lu) e{uj) , 
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Fig. 10. Stress relaxation analysis for the SiH material. Graph 
(a) shows the DMA Master Curve, from Fig. ^(a), plotted 
against the inverse frequency The region of long times 

is expanded in Graph (b), where the real relaxation data from 
Fig. pl(a) are shifted along the time-axis to superpose on the 
same Master Curve. The shift factors required for this super- 
position are shown as square symbols in Fig. |^(b). 
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Fig. 11. Stress relaxation analysis for the SiMC material. 
Graph (a) shows the DMA Master Curve, from Fig. ^(a), plot- 
ted against the inverse frequency uj~^ . The region of long times 
is expanded in Graph (b), where the real relaxation data from 
Fig. H(b) are shifted to superpose on the same Master Curve. 
The shift factors required for this superposition are shown as 
square symbols in Fig. rab). 



and assuming the linear response function to decay as a 
power law, at long times G{t) « Goi~°, one can easily 
calculate the corresponding estimate of the complex shear 
modulus: 



G 



(3) 



Apart from a constant coefficient, this is directly the in- 
verse function of G{t). 

Let us now take the real-time stress relaxation data 
sets and attempt the time-temperature superposition on 
them, building a Master Curve by shifting along the time- 
axis. This operation works for all temperatures, but we 
only present three data sets in each of Figs. |l^(b) and 
|ri|(b), to allow a clear view of the underlying line, which 
is the DMA Master Curve for inverse frequencies. The 
matching is quite remarkable, especially if we take a look 
at the shift factors (now for the time-shifts, and so with the 
opposite sign!) required to bring the relaxation data onto 
this Master Curve. Figures ^(b) and ||(b) show the fre- 
quency shift factors and their WLF analysis, but also the 
square symbols that represent the (inverse) time shift fac- 
tors for the static relaxation data. (Recall that the Master 
Curves were built for the reference temperature of around 
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20°C). The agreement and the trend in these respective 
shift factors is not bad. 



5 Discussion 

First of all, we must recognise the difference of our results 
from some literature data. Gallani et al. and Weilepp et 
al. 10,^ were both pointing out that the nematic-isotropic 
transition is trivial for the dynamic mechanical measure- 
ments. The reason could be that polydomain nematic elas- 
tomers were used, although we have earlier seen the effect 
of dynamic soft elasticity in a variety of polydomain elas- 
tomers . The drop in G' below Tni is smaller in polydo- 
main samples, but it only disappears in the monodomain 

confirms. 
, reported 



V-geometry, as in fact the more recent studyjll | 
It is also interesting that many authors, e.g. 
their success in time-temperature superposition across the 
nematic-isotropic phase transition, which of course is re- 
lated to the absence of pronounced soft drops in G' values. 
We can only suppose that a different synthetic procedure 
of their materials led to the very non-soft (highly semi- 
soft) nematic networks. 

Returning to our own results, the two main results of 
this work are the building of Master Curves and iden- 
tifying the temperature regimes where the superposition 
follows the WLF and where it fails altogether, and es- 
tablishing the direct analogy between the time and the 
frequency domains. In the first part, the most important 
observation is the substantial difference in rubber-plateau 
(low-frequency) moduli between the nematic and isotropic 
phase, in the shear geometry when the director rotation 
is induced. This effect of dynamic soft elasticity can lead 
to a number of technological advances, e.g., in the area of 
mechanical vibration damping p6[ (exploiting the record- 
high values of the loss factor over a broad range of temper- 
atures and frequencies), or even suggesting a novel field of 
polarised acoustics (using the fact that only certain, soft, 
shear- wave geometries are strongly attenuated). 

Comparing the results of static stress relaxation with 
those of dynamic mechanical measurements is intriguing. 
Regarding the ability to quantitatively predict the static 
relaxation behaviour from the independently (and rela- 
tively quickly) acquired DMA data, or vice versa, by in- 
verting the Master Curve at a specific reference tempera- 
ture, we have to remain cautious. Although the possibility 
is highly attractive and the procedure has worked spectac- 
ularly well in our case, one needs to verify it on different 
materials, including a check of how the classical isotropic 
rubber behaves in this respect. 

In our case, the Master Curve inversion has been suc- 
cessful and predicts that, for instance, a "traditional" 
|20|| siloxane side-chain nematic elastomer SiH, crosslinked 
in a way that allows the soft elasticity, would be able to re- 
lax under a constant strain and reach its mechanical equi- 
librium (at least approximately) after some lO^s at room 
temperature, with a very little residual relaxation after 
that. This is somewhat different from the results on 
a very similar, but polydomain nematic elastomer, where 



the relaxation have not reached its end even after extrap- 
olation to astronomical times. Analogous continuing re- 
laxation was reported in other, chemically very different 
polydomain nematic rubbers ||l^. The effect of internal 
mechanical constraints imposed by the domain interfaces 
is clearly dominating long-time regime, but is absent in our 
present monodomain network. In contrast, the composite 
nematic elastomer SiMC, albeit monodomain, also shows 
a non-ending relaxation. We assume, by analogy with the 
previous work on relaxation in non-crosslinked MC poly- 
mer that internal constraints on "free" stress relax- 
ation are now imposed by the hairpin folds of long nematic 
main-chain polymer strands. 

In summary - dynamic-mechanical analysis have been 
performed using two different monodomain nematic elas- 
tomers having different degrees of anisotropy and internal 
microstructure. Master Curves have been built between 
the glassy state and the nematic-isotropic phase transi- 
tion, where the modulus saturates on the low-level soft 
plateau. Above T„i the modulus rises substantially, since 
the internal relaxation is no longer able to reduce the elas- 
tic response - and the further time-temperature superpo- 
sition fails. The dynamics of these elastomers is dominated 
by power laws (non-universal, depending on the material 
and temperature), which is confirmed by the successful op- 
eration of the Master Curve inversion to describe the static 
stress relaxation. Power laws (in time) of stress relaxation 
match very well with the corresponding frequency depen- 
dence of the dynamic modulus in the appropriate range of 
temperatures. 
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